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Abstract 

In the formalism of measurement based quantum computation we start with a given 
fixed entangled state of many qubits and perform computation by applying a sequence 
of measurements to designated qubits in designated bases. The choice of basis for later 
measurements may depend on earlier measurement outcomes and the final result of the 
computation is determined from the classical data of all the measurement outcomes. 
This is in contrast to the more familiar gate array model in which computational steps 
are unitary operations, developing a large entangled state prior to some final measure- 
ments for the output. Two principal schemes of measurement based computation are 
teleportation quantum computation (TQC) and the so-called cluster model or one-way 
quantum computer (1WQC). We will describe these schemes and show how they are 
able to perform universal quantum computation. We will outline various possible rela- 
tionships between the models which serve to clarify their workings. We will also discuss 
possible novel computational benefits of the measurement based models compared to 
the gate array model, especially issues of parallelisability of algorithms. 



1 Introduction 

Many of the most popular models of quantum computation are direct quantum generali- 
sations of well known classical constructs. This includes quantum turing machines, gate 
arrays and walks. These models use unitary evolution as the basic mechanism of infor- 
mation processing and only at the end do we make measurements, converting quantum 
information into classical information in order to read out classical answers. In contrast to 
unitary evolution, measurements are irreversibly destructive, involving much loss of poten- 
tial information about a quantum state's identity. Thus it is interesting, and at first sight 
surprising, that we can perform universal quantum computation using only measurements 
as computational steps EJ El Ej • These "measurement based" models are especially 
interesting for fundamental issues: they have no evident classical analogues and they of- 
fer a new perspective on the role of entanglement in quantum computation. They may 
also be interesting for experimental considerations, suggesting a different kind of computer 
architecture and offering interesting possibilities for further issues such as fault tolerance 



We will discuss two measurement based models. Firstly we'll consider teleportation 
quantum computation (TQC). This is based on the idea of Gottesman and Chuang |2] of 
teleporting quantum gates, and was developed into a computational model by Nielsen, Le- 
ung and others [HI E] ■ Secondly we will consider the so-called "one way quantum computer" 
(1WQC) or "cluster state computation" of Raussendorf and Briegel jSJE]. Then we will 
discuss various possible relationships between these two models and finally consider possible 



computational benefits of the measurement based formalism as compared to the quantum 
gate array model. 

The following notations will be frequently used. The Pauli operators are 

*-(;;) *-(;.;) j 

© will denote addition modulo 2. The controlled-NOT gate CX is defined by 

CX\i)\j) = \i)\i®j) i,i = 0, 1. 
The controlled phase gate CZ is defined by 

CZ\i) \j) = (-lf\i) 

In contrast to the CX gate, CZ is symmetrical in the two input qubits. The Hadamard 
operator is 

and the Hadamard basis states are 

l±) = -L(|o>±|i». 



The Bell basis states are 



|B 00 ) = TaGOO) + |11» = I®I\B 
\B 01 ) = ^(\01) + \10))=X®I\B 
\B W ) = ^(|00) - |11)) = Z®/|5, 
|S al ) = ^(|01)-|10))=iF® I|£ 00 }. 
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Succinctly we have 



\B cd ) = Z c X d ®I\B l 



oo 



We will also use the maximally entangled state that combines the Hadamard and standard 
bases in its Schmidt form: 

|ff> = (| ) |+) + |1) |-)) = ± (|+) |0) + |-) |1» = CZ\+) |+) . 
2 Teleportation based quantum computing 

Recall standard teleportation [J as depicted in figure 1 and explained in the caption. 

Now introduce the "rotated Bell basis" denoted B(U) = {\B(U)cd}} where for any 1- 
qubit gate U : 

\B(U) cd ) = U^I\B cd ). 

A simple calculation shows that if we perform such a rotated Bell measurement (instead of 
the standard one in figure 1) then the teleported state at 3 is X d Z c U \ip) i.e. the gate U 
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Figure 1: Standard quantum teleportation. The labels 1,2,3 represent three qubits. 1 is in state 
| a) and 23 are in state |£>oo)- A Bell measurement is performed on 12. If the outcome is \B c d) then 
qubit 3 acquires the state X d Z c \a). 





Figure 2: Projecting the three qubit state \a) 1 \4>) 2 3 onto the maximally entangled state \<fr) 12 results 
in a single qubit state at 3 given by 4 |a) 3 (see lemma 1). 



has been applied to \tp) via this teleportation 2 j. A particularly neat and general way of 
performing this calculation is the following. For any dimension d consider any maximally 
entangled state, written in its Schmidt form as 

d-l 

Consider the mathematical projection depicted in figure 2. 

Lemma 1: The projection of \a) 1 |0) 23 onto \4>} 12 results in the state -? \a} 3 at qubit 3. 
Proof: Let \a) = Y2 a j \ Then the projection is 

i jk ijk k 

Next introduce the "rotated | <j>) state" : for any unitary operator U in d dimensions define 

\<f>(U)) = tf®I\<f>). (1) 

Lemma 2: The projection of |a) 1 |<^) 23 onto |0(C/)) 12 results in the state -ill |a) 3 at qubit 
3. 

Proof: For any states \a), \b) if we write It^'ct) = \a) then <^J7^a|6) = (a\Ub). Hence 
lemma 1 with \(j)} 12 replaced by \4>(U)) 12 immediately gives the result. □ 
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Figure 3: Application of CZ by teleportation with 8 dimensional measurements. The wiggly 
line connecting 12 denotes an input 2-qubit state \ip). The lines labelled H denote the maximally 
entangled state \H). Introduce the 3-qubit Bell measurement corresponding to the basis { |000) ± 
| 111) , |001) ± |110) , |010) ± 1 101) , |100) ± |011)}. If this measurement is performed on qubits 135 
and 246 then it can be verified that qubits 78 acquire the state (Pi ® Pj){H <£> H)CZ \ ip) where the 
Pauli operators Pi and Pj depend on the measurement outcomes. 



Now the basic idea is to regard the 12-projections in lemmas 1 and 2 as being outcomes 
of a projective measurement applied to qubits 1 and 2 i.e. we wish to choose a set of d 2 
unitaries U{ such that {\4>(Ui))} form an orthonormal basis of the 12 space. For qubits 
(d = 2) the set {/, X, Z, XZ} gives standard teleportation and the set {IU, XU, ZU, XZU} 
gives the rotated Bell basis, reproducing our previous procedure of "gate teleportation" to 
construct X d Z c U at qubit 3. 

Note that lemmas 1 and 2 apply in general dimension d so we can apply 2-qubit gates 
such as CZ via measurements by teleportation in dimension d = 4. For example using 
\4>) = | i?oo) I -Boo) an d the 16 operators C/y = (Pi ® Pj)CZ where Pi and Pj range over the 
four standard Pauli operators, we can check that {U}j ®I\<p}} is an orthonormal set and 
the output teleported state is (Pi <g> Pj)CZ for any input 2-qubit state \ip). 

This method of applying CZ requires a 16 dimensional Bell measurement but with more 
subtle means we can achieve the result with smaller measurements. An illustrative example 
(taken from |18j ) is shown in figure 3. 

Remark: Returning to \4>(U)) in eq. (^Q), a set of d 2 operators {Ui} will have the 
corresponding set {\4>(Ui))} being orthonormal iff Tr(UiUj) = 5ij i.e. {Ui} is a so-called 
unitary operator basis. For general dimension d many such sets exist, each corresponding to 
a teleportation scheme. Werner |Sj has described a method for constructing a large number 
of explicit inequivalent examples. However we can go even further and choose a set of n > d 2 
operators Ui such that the set {ki \4>(Ui)) (<f>{Ui)\} (for some chosen constants ki > 0) form 
the elements of a (rank 1) positive operator valued measure (POVM) i.e. we ask that 



The case of n = d 2 (with all ki's then necessarily equal) reproduces projective measurements 
as above but for n > d 2 we obtain fully valid teleportation schemes in which the Bell 
measurement is replaced by a POVM. Even in the qubit case d = 2 many examples exist 
and in fact we can have unboundedly large n. (We will not digress here to discuss explicit 
examples of such constructions. )□ 



n 



Y J h\^{u i ))((l>(u i )\=i d ®i d . 
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Thus in summary so far, if we have a pool of maximally entangled states we can apply any 
unitary gate U to any (multi-qubit) input state \ip) by measurements alone. A significant 
annoyance (and see more about this later) is that we do not get the exact desired result 
U but instead get PU \ where P is some Pauli operation (on each qubit) depending on 
the measurement outcome. This is the residue of the randomness of quantum measurement 
outcomes in our computational formalism. 

To perform universal computation it suffices to be able to apply gates from any conve- 
nient universal set. We introduce notation for further gates, x and z rotations are defined 
by 

R x (0) = e- wx R z (e) = e- idz . 
The following gate will be significant for 1WQC later: 

W V) = 7f ( \ S ) = Hm where m = ( o J?) • (2) 

Any 1-qubit gate U (up to an overall phase) can be decomposed as (Euler angles): 

U = R x (C)RMRx{0 for some £,t],C 

and also ^Hl as 

U = W(0)W(e 1 )W(9 2 )W(9 3 ) for some 0i,0 2 ,03 • 
It is known that CZ together with all 1-qubit operations is a universal set so the sets 
{CZ,R x (e),R z ((f)) all 9,<t>} and {CZ,W{6) all 6} are also universal. 



3 Adaptive measurements 

We wish to perform arbitrary sequences of gates from a universal set by measurements. For 
simplicity consider 1-qubit gates: to perform . . . U3U2U1 \ip) we successively teleport the 
three gates but instead of the desired result we get . . . P3U3P2U2P1U1 \ip) where Pi, P2, P3 . . . 
are Pauli operators depending on the measurement outcomes. To deal with this awkward- 
ness we introduce a new feature: adaptive choices of measurements. Let us assume that 
all 1-qubit operations are x or z rotations. The following Pauli "propagation" relations are 
easily verified: 

R X (6)X = XR X (6) R X {6)Z = ZR x {-6) 
R Z (9)X = XR z (-9) R Z {6)Z = ZR Z {6). 

(Alternatively we could use only W / (0)'s with relations W(0)X = ZW(-6) and W{0)Z = 
XW{6).) Suppose we want to do . . . R z {(3)R x {a) \tp). The first Bell measurement with out- 
comes a, b gives Z a X b R x (a) \ip). Now commuting Z a X b to the left through R Z (P) causes (3 
to change to (—\) h [3. So if instead of R z {(3) we next applied the measurement for R z {{—l) b f3) 
instead, then 

R z ((-l) b P)Z a X b R x (a) |V>> = Z a X b R z ((3)R x (a) |V>> 

as wanted. This second measurement with basis adapted to measurement outcome b has an 
outcome c, d say, and the state is now Z a+c X b+d R z (f3)R x (a) \ ip). Continuing in this way, 
adapting measurement bases to earlier measurement results, we get 

. . . Xjp*2% a X? 1 2% 1 (the correct wanted U) k/>) . 
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Here is generally a state of many qubits and Xi, Zj are Pauli operations on the i qubit. 
The indices mi, rii are accumulations (actually bit sums mod 2) of measurement outcomes. 

The same idea applies to the 2-qubit gate CZ too where the situation is even better. 
We have the propagation relations: 

CZ{Z®I) = (Z®I)CZ CZ{X®I) = {X®Z)CZ 

(and similarly for X and Z acting on the second qubit on the LHS's, as CZ is symmetrical). 
Thus we can propagate Pauli operators through CZ while keeping CZ the same i.e. no 
basis adaption is required! 

Remark: More generally for any 1-qubit U we can achieve Pauli propagation by 

UX = X(XUX) UZ = Z(ZUZ) 

where we adapt U to change into XV X or ZUZ. But the latter operators look very 
different from U and we have (unnecessarily) also preserved the actual identity of the Pauli 
operations in this propagation. The actual relations we used above for CZ, R x (6), R z {6) 
and W{6) exploit the opportunity of allowing the Pauli operations to change while keeping 
the adapted operation similar to the original one (e.g. differing only by a sign in the angle). 
□ 

If U is the total unitary effect of a gate array on a multi-qubit state \ip) then using the 
above methods we are able to generate a state of the form . . . X™ 2 Z^X™ 1 Z™ 1 !/ \tp). In 
a quantum computation we finally measure (some qubits of) U in Zj-bases {|0) , |1)} 
and the presence of the Pauli operations X^Z™ 1 cause no problems: Z™° has no effect on 
the measurement outcomes and X™ % requires only a simple reinterpretation of the output 
results - for a single qubit, if U\ip) = a |0) + b |1) then XU = a |1) + b |0). Hence we 
simply need to reinterpret measurement outcome ki (0 or 1) as h L © mj where mj is the 
corresponding X Pauli exponent. 

We note a rather curious feature here: these "final" Z measurements are never adaptive 
(being fixed as Z measurements) so they can always be performed first before any of the 
other measurements have been implemented! i.e. the output of the computation can be 
measured before any of the computation itself has been conducted and the subsequent 
measurement outcomes simply serve to alter the interpretation of those Z-measurement 
outcomes! 

4 Parallelisable computations: Clifford operations 

In our measurement based models, distinct measurements always apply to disjoint sets of 
qubits. Hence they all commute as quantum operations and if the measurement basis choices 
were not adaptive then we could do all the measurements simultaneously, in parallel. The 
necessity of adaptive choices arose from the Pauli propagation relations, but some operations 
are special in this regard - the so-called Clifford operations on n-qubits. 

Let us introduce the Pauli group V n on n qubits, defined as the group generated (under 
multiplication) by n-fold tensor products of ±1, ±il, X and Z. For example Vz has elements 
such as X % Z % Z, —iZ %Y %I, 1 1 X etc. An n-qubit unitary operation C is defined 
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to be a Clifford operation if CV n C^ = V n i.e. CV n = V n C i.e. for every Pauli operation 
P G V n there is another P' £ V n such that CP = P'C. Hence for any Clifford operation 
C we can propagate Pauli operations across C while C stays the same i.e. no adaption 
is needed (but the Pauli operation generally changes). For each n the Clifford operations 
evidently form a group, called the Clifford group. 

We have already noted that CZ is a 2-qubit Clifford operation. Also the Hadamard 
operation has the Clifford property since HX = ZH and HZ = XH. (For any prospective 
C we need only check the propagation of X and Z at each qubit to verify the full Clifford 
property.) Indeed we can give an explicit description of the Clifford group on n qubits 
Introduce the 7r/4 phase gate: 



Then we have: 

Theorem: 9 The Clifford group on n qubits is generated by Z, H, P^u and CX acting 
in all combinations on any of the qubits (i.e. arbitrary arrays of these gates). □ 

Hence any array of Clifford gates Ck . . . C2C1 \ip) (\ip) of n qubits) may be imple- 
mented in TQC in one parallel layer of measurements. We get a result of the form 
PfcCfc . . . P2C2P1C1 where P{ are all Pauli operations on the qubits. Commuting them 
all out we get X^ n Z^ n . . . JT" 1 Z* 1 (Cfc . ..Ci)\tp) where the indices a;, bi depend on the mea- 
surement outcomes and the Clifford propagation relations. The collection of maximally 
entangled states used in all the teleportations can also be manufactured in parallel (e.g. 
apply CX's to many pairs |0) |0)) so the entire quantum process requires only a constant 
amount of quantum parallel time for any n, in contrast to the corresponding gate array 
whose depth generally increases with n. 

However there is a further subtle point here: in addition to the constant parallel time 
of the quantum process, the computation of the Pauli exponents Oj, bi requires a further 
classical computation, which is actually the bitwise sum of selected measurement outcomes. 
Hence this classical computation can be done as a parallel computation of log depth: to 
sum k bits i\ . . . %\~ we first sum all pairs i\ © 12, 13 © 24, • • • in parallel and then pairs of the 
results etc. At each stage the number of bits is halved so we require log k layers to reach the 
final result. Returning to the gate array of Clifford operations, it is known that any Clifford 
operation on n qubits can be represented as an array of O(logn) depth i.e. the same as the 
total quantum plus classical depth of the full measurement based implementation. But the 
virtue of the measurement based approach is to take a fully quantum process (in this case, 
the Clifford gate array of log depth) and recast it as a quantum process of constant depth 
plus a classical computation (of log depth in our case) i.e. we separate the original process 
into a quantum and classical part while suitably "minimising" the quantum part. This kind 
of restructuring is significant in a scenario where quantum and classical computation are 
regarded as separate (or even incomparable) resources (c.f. further discussion in section 
|H] below) and we can ask interesting questions such as: what is the "least" amount of 
quantum "assistance" needed to supplement classical polynomial time computation in order 
to capture the full power of quantum polynomial time computations? 
Remark: It is known that arrays of quantum Clifford operations can be classically effi- 
ciently simulated. This is the Knill-Gottesman (KG) theorem jlUj which asserts the follow- 
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Figure 4: Construction of the cluster state. The dots represent a grid of |+) states and a line 
connecting a pair of dots represents the application of CZ to the corresponding pair of qubits. 



ing: consider any array of Clifford gates on n qubits, each initialised in state |0). Let V 
be the probability distribution resulting from measuring (some of) the output qubits in the 
Z basis. Then there is a classical (probabilistic) computational process which runs in time 
polynomial in n, which also has V as its output distribution. Furthermore according to 
jllj this classical simulation of any Clifford array's output can always be performed in log 
depth. In view of this, one may question the significance of our result above, that any such 
Clifford array can be reproduced as a classical (log depth) process plus a quantum process 
of constant depth. However there is an essential difference in the two representations of the 
Clifford array: the KG simulation results in a purely classical output (sample of V) whereas 
our classical-quantum separated simulation results in a quantum state as output i.e. we are 
simulating the quantum process itself rather than just the classical output of some measure- 
ment results. Indeed our result makes an interesting statement about quantum properties 
of Clifford arrays, viz. that the essential quantum content can be "compressed down" into 
constant quantum depth, which is not provided by the KG result. □ 

5 The "one-way" quantum computer 

We now move on to describe our second measurement based computational model - the 
1WQC of Raussendorf and Briegel [SJ At first sight it looks rather different from TQC 
but we will see later that the models are in fact very closely related. 

Consider a rectangular (2 dimensional) grid of |+) states as shown in figure^]. We apply 
CZ to each nearest neighbour pair (in horizontal and vertical directions). These CZ's all 
commute so for any grid size they can all be applied in parallel, as a process of constant 
quantum depth. The resulting state is an entangled state of many qubits, called a cluster 
state. We will also use one dimensional cluster states constructed in the same way, but 
starting from a one dimensional array of |+) states. 

The 1WQC is based on the following facts that will be elaborated below. Any quantum 
gate array can be implemented as a pattern of 1-qubit measurements on a (suitably large two 
dimensional) cluster state. The only measurements used are in the bases M z = {|0) , |1)} 
and M{9) = {|0) ± e t8 |1)} for some 8.9 = corresponds to an X basis measurement. 
Measurement outcomes are always labelled or 1 and they are always uniformly random. 
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Figure 5: 1WQC implementation of a 1-qubit unitary U on \ip). £, rj, and ( are Euler angles for U. 
The leftmost qubit, denoted by a star, is set in state l^) and extended by a row of four |+) states 
denoted by dots. CZ operations are then applied, denoted by connecting lines. Next, measurements 
are applied in the designated bases with outcomes Sj. Hence the measurements must be carried out 
adaptively from left to right. As a result of this process the rightmost (unmeasured) qubit is left in 
state X S2+Si Z Sl+S3 U\ip). 



As in teleportation, quantum gates are implemented only up to Pauli corrections X a Z b 
where a and b depend on the measurement outcomes. (These Pauli corrections are called 
bi-product operators in the 1WQC literature). Hence we'll get the same feature of adaptive 
measurements that we saw in TQC. The name "one-way" quantum computer arises from 
the feature that the initial resource of the pure cluster state is irreversibly degraded as the 
computation proceeds in its layers of measurements. 

To illustrate these ingredients we give some explicit examples of measurement patterns 
for 1-qubit gates (where 1-dimensional cluster states suffice). Our first example is taken 
from f(3j. We noted previously that any 1-qubit U can be expressed as 

U = R x (C)RMRx(0 for some ^,V,C- 

To apply U to \ r ip) by the 1WQC method we start with \ r ip) in a line with four |+) states, 
as shown in figure [5J (Later we will see how to eliminate explicit use of here, starting 
with only |+) states). Entangle all neighbouring pairs with CZ and subsequently measure 
qubits 1,2,3 and 4 adaptively in the bases shown in figure |SJ Then it may be shown that 
qubit 5 acquires the state X S2+S4 Z Sl+S3 U \ where Sj is the outcome of the measurement 
on qubit i. Note that the measurement bases are adaptive and in the above pattern the 
measurements must be performed in numerical sequence. 

As a second example consider the pattern of figure Then qubit 3 acquires the state 
Z Sl X S2 so we have a process which is very similar to teleportation from 1 to 3. The two 
X measurements may be done in parallel but we can also consider this process as a sequence 
of two identical steps: given a state adjoin |+), entangle with CZ and then X-measure 
the first qubit, giving an outcome s. Qubit 2 is then left in state X s H \ tp). Thus the two- 
step chain of figure El can be analysed as (X S2 H)(X Sl H) \if)) = X S2 Z Sl HH \tp) = X S2 Z Sl 
(where we have used the Pauli propagation relations for H). 

Figure [7| shows a single step operation for the general measurement basis M{6). Qubit 
2 is then left in state X Sl W(— 6) \tp) (with W(8) as defined in eq. (J2J)). Indeed if we are not 
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Figure 6: See caption of figure [3] for explanation of pictorial notations. After the two X measure- 
ments, qubit 3 is left in state Z Sl X S2 \ip) where \tp) was the input state at qubit 1. 
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Figure 7: Effect of a single M{6) measurement in 1WQC. Qubit 2 is left in state X Sl W(-6) \ip). 



concerned with issues of parallelisability then any 1-dimensional measurement pattern may 
be viewed as a sequential application of the process of figure [7| applied repeatedly. This 
shows that the unitary operation W(0) plays a fundamental role in 1WQC. For example 
in figure [3 we can reorder all operations as follows: (entangle 12, measure 1), (entangle 23, 
measure 2), etc. This is because the 12-entangling operation and first measurement both 
commute with all subsequent entangling operations and measurements. 

In a similar way it is now straightforward to see how gates may be concatenated. Suppose 
we wish to apply U\ and then U2 to an input qubit \ifi). Measurement pattern 1 for U\ has an 
output qubit (e.g. qubit 5 in figure^ which is the input qubit for the measurement pattern 
of 1/2- But all measurements in pattern 1 commute with all entangling operations of pattern 
2. Hence we can apply all entangling operations (for both patterns) first, to get a longer 
single cluster state and then apply the measurements. Furthermore some measurements in 
pattern 2 could even be performed before those in pattern 1 if their basis choice does not 
depend on pattern 1 outcomes. 

Finally we can eliminate the input state from the above descriptions, to get a 
formalism based entirely on a starting state that's a fully standard cluster state (of slightly 
longer length): since we can implement any 1-qubit U we can take our input starting state 
to be |+) and prefix the desired process with an initial measurement pattern for a unitary 
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Figure 8: Application of CZ gate in 1WQC. The 2-qubit input state l^in) is placed at sites labelled 
ini and in2. Dots denote |+) states and connecting lines denote application of CZ for cluster state 
generation. If the measurement pattern shown at the right is applied at the sites, then only sites 
outi and out2 remain unmeasured and contain (Pi ®P 2 )CZ \^pi n ) where P\®P2 is a Pauli operation 
that depends on the measurement outcomes. 
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Figure 9: Z measurements are used to eliminate unwanted sites such as A. 



operation that takes |+) to 

Above we have considered only 1-qubit gates but the formalism may be generalised 
(using 2-dimensional cluster grids) to incorporate 2-qubit gates. For universal computation 
it suffices to be able to implement just the CZ gate (in addition to 1-qubit gates). An 
explicit measurement pattern for CZ is shown in figure El It should also be noted that 
measurement patterns are not unique and subject to various approaches for their invention. 



5.1 Role of Z measurements 

The CZ gate requires use of a 2-dimensional grid whereas 1-qubit gates require only 1- 
dimensional clusters. Hence measurement patterns for general gate arrays, implemented on 
a suitably large 2-dimensional cluster, will generally have some extraneous sites not used in 
the measurement patterns. Z-measurements are used to delete such extraneous sites. 

To illustrate the principle consider a cluster state with irregular shape in figure El 
Suppose we wish to use only the linearly arranged sites 1,2,3 i.e. we want to delete site A. 
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Consider Za, a Z measurement at site A, with outcome kA- To see its effect recall that 
the pattern in figure is obtained by applying CZ operations to a set of four |+) states, 
located at the sites. But measurement Za commutes with CZyi and CZ23 so starting 
with the four |+) states we can first do CZa2 and then Za before CZ\2 and CZ23. But 
CZa2 |+) |+) = |0) |+) + |1) |— ) so after the Za measurement (with outcome fe^) the sites 
1,2,3 contain |+) |(— l) fc - 4 ) |+) i.e. the effect of the Za measurement is to have used the 
Hadamard basis state |(— l) fc - 4 ) at site 2 instead of the standard |+) (with no site at A) 
and then entangle as usual. Next consider the basic 1WQC single step of figure [7] with |— ) 
instead of |+) at site 2: we start with |— ), apply CZ, then do M{8) at qubit 1 for output 
at qubit 2. But this is the same as the following: start with |+), then (i) apply CZ, (ii) 
do M(6) on qubit 1, then (iii) perform the unitary operation I\®Zi- This is identical to the 
previous process because (iii) commutes with (i) and (ii) so we can change the order to (iii), 
then (i) then (ii) and note that Z2 |+) = |— ). A similar argument applies to each neighbour 
of A if there is more than one. In summary, we see that a Za measurement (outcome kA) 
at a site A adds in an extra Z kA Pauli correction at all neighbouring sites of A, in addition 
to the usual Pauli operations arising from measurement patterns in a standard cluster state 
that had site A absent from the start. 

Z measurements also have a second role: as in TQC they are applied to a final state of a 
computation to produce the classical output results. So just as in TQC we have the curious 
feature that these "final" Z measurements are always non-adaptive and can be applied first, 
before any of the computation itself has been implemented! 

To briefly summarise the 1WQC model, we have seen that any quantum gate array 
can be translated into a pattern of 1-qubit measurements on a suitably large 2-dimensional 
cluster state. Choices of measurement bases are generally adaptive i.e. possibly depending 
on previous measurement outcomes. Thus the 1-qubit measurements are organised into 
layers and the measurements within each layer can be done simultaneously in parallel. The 
output qubits (always measured in the Z basis) can always be done in the first layer. This 
temporal sequencing "cuts across" the temporal sequence of the original gate array - all 
gates are generally done "partially" in each layer and simultaneously built up as the layers 
accumulate. Just as in TQC, arrays of Clifford operations can always be fully implemented 
with only one layer of measurements. 

6 Further features of 1WQC 
6.1 A further parallelisability result 

Any polynomial sized quantum gate array can be implemented in 1WQC using at most a 
polynomial number of measurement layers (c.f. section |H] later). Also we have seen that any 
array of Clifford operations can be implemented with just one layer of measurements. This 
suggests the following interesting question: which classes of quantum gate arrays can be 
implemented in 1WQC with constraints on the number of measurement layers e.g. using 2 
or 3 or a logarithmic number of layers? Does the latter include all polynomial time quantum 
computation? Although very little is known about such questions, we have the following 
result of Raussendorf and Briegel E]- Its proof depends on more subtle properties of 
Pauli propagation relations for particular operations. 
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Figure 10: The process above leaves qubit 3 in state X S2 Z Sl R x ((— l) Sl 8) \^p). Thus to implement 
a given x rotation the sign of the angle 9 is adaptive, depending on the measurement outcome s±. 



Theorem: Any gate array using gates from the set {CX, R x (6) all 9} or from the set 
{CX, R z (9) all 9} can be implemented with just two measurement layers. 
Remark: Neither of these sets is believed to be universal although it is known that CX 
with all y-rotations is universal jl.Sj . 

Proof of theorem: We give a proof for x-rotations. (The case of z-rotations is similar). 
We use the following three (easily verified) facts: (i) any single M(9) measurement (e.g. 
as shown in figure [7J) generates only I or X Pauli corrections. They may become Z's only 
after propagation through other gates, (ii) when Pauli operations are commuted across the 
Clifford operation CX, X propagates only to X's (and Z propagates only to Z's) although 
the Pauli operator may spread from one qubit onto two qubits. (iii) The R x gate may be 
implemented in 1WQC using the measurement pattern shown in figure ITU1 

Now consider any gate array . . . G3G2G1 where each Gi is a CX or an x-rotation gate. 
Set up its corresponding measurement pattern (using figure ITU1 for each x-rotation gate). 
CX is a Clifford operation and involves no adaptive measurements. In the first measurement 
layer we perform all CX pattern measurements and all X measurements of the R x gates (i.e. 
all qubit 1 measurements in figure ITU)) . This produces some extraneous Pauli operations and 
leaves only the M(—9) nodes unmeasured (i.e. all qubit 2's in figure IT0|) . Next commute 
all these Pauli's out to the left hand end of the gate array. This commutation leaves all 
CX's unchanged (as CX is Clifford) but some x-rotation angles acquire an unwanted minus 
sign (c.f. propagation relations for x-rotations given previously). Thus reset these altered 
angle signs so that the M(—9) measurements again give the correct designated gates i.e. 
make the adaptive choice of bases for the next measurement layer. Finally in the second 
layer, perform all the M{—9) measurements. We again get some further extraneous Pauli 
operations generated. By (i) we get only / or J, but these can be harmlessly commuted 
out to the left since X commutes with x-rotations and CX preserves X's in its propagation 
relations (i.e. no Z , s are generated). □ 

6.2 Non-universality of one-dimensional clusters 

We have seen that 1WQC with two-dimensional cluster states is universal for quantum 
computation. Nielsen and Doherty ^1] have shown that any 1WQC process based on only 
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Figure 11: A generalised 1-dimensional cluster state \$ n }- Starting with a row of n qubits each in 
state |0) we sequentially apply L/j^+i to qubits i,i + 1 as shown in the above gate array having a 
ladder-like structure. 



1-dimensional cluster states can be simulated classically efficiently i.e. in polynomial time in 
the number of qubits. Thus the universality of any such model would imply that quantum 
computation is no more powerful than classical computation. 

The argument may be paraphrased as follows (ignoring technical issues of precision of 
the simulation). Consider any state of the form shown in figure ITT1 (The standard cluster 
state is obtained by choosing each C^+i to be CZ(H®H)). Consider now any sequence of 1- 
qubit measurements such that for later measurements, the choice of measurement basis and 
even the choice of qubit used, may both depend on outcomes of earlier measurements. Then 
the whole process may be classically simulated in polynomial (in n) time i.e. the resulting 
probability distribution of outcomes may be sampled by classical means in polynomial time. 

The proof runs as follows. At a general stage suppose measurements on lines a,b, . . . ,p,q 
have already been simulated and the chosen sample outcomes are a' ,b' , . . . ,p' , q' respectively. 
Even though these measurements may have been chosen adaptively, once the outcomes have 
been specified (as a' , 6', . • • iP q') we can perform the measurements in any order we wish 
to compute the joint probability P(a' , 6', . • • )P '> <f) of the designated outcome string. Now 
suppose the next measurement is on line k. To sample its outcome distribution we need to 
know 

Here a' ,b' , . . . ,p' ,q' have their fixed values and we take k' = 0, 1 separately. Having com- 
puted this probability we sample k! to give a definite value, and continue in the same way 
for the next measurement on some line, I say. 

To show that this whole process is efficient, we need only show that P(a! , b' , . . . ,p', q') 
can be computed in poly(n) time for any chosen set of values a', b\. . . ,p' , q' , on any chosen 
set of lines. Without loss of generality suppose these lines are listed in increasing order of 
occurrence from line 1 at the top. Note that we can regard a measurement on any line, a 
say, as occurring immediately after U a ^ a +i and before C/ a+l a+2 . Let p m denote the reduced 
state of any line m at the position in between the application of J7 TO -i,m and U m . m+ \, and 
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let cr m) m+i denote the application of U m>m +i to p m (8> |0) (0|. 

Starting at the top we compute a 12 as U12 on |0) |0). Then compute P2 by partial trace, 
tracing out system 1 from oyz- Then compute 023 as XJ2Z on p2® |0) (0| etc. continuing until 
o a ,a+i where the first measurement occurs. At this point we do not compute p a+ \ as the 
partial trace above but instead, apply to a a ,a+i the projector corresponding to measurement 
outcome a', obtaining a subnormalised state p a +\ on line a + 1. In fact Trp a+ i is the 
probability of getting outcome a 1 in the a- line measurement. We continue in this way 
computing the reduced states of the lines successively (using the measurement projector at 
any measured line and partial trace at any unmeasured line) until we have applied the last 
(q') measurement projector. The trace of this final resulting state is then P{a! ,b' , . . . ,p' , q'). 
At each stage of this calculation we need to hold the state of at most two lines (i.e. not 
n lines with its exponentially large description) and we pass through the ladder n times, 
once for each successive measurement. Hence the whole calculation is completed in time 
polynomial in n. This feature of the calculation is a consequence of the special ladder-like 
structure of figure ^2 and it does not apply to general gate arrays. 

7 Relationships between the TQC and 1WQC models 

The two models have several similarities - both are based on measurements as compu- 
tational steps and both have the awkward feature of supplementing desired gates with 
unwanted Pauli operations. But there are also some essential differences: TQC uses (Bell) 
measurements on 2 or more qubits whereas 1WQC uses only 1-qubit measurements. 1WQC 
starts with a cluster state having multi-partite entanglement across all the qubits whereas 
TQC is based on a state comprising only bipartite entangled pairs. 

The models can be related in several different ways. We will discuss three of them. The 
relationships serve to improve our understanding of 1WQC and its measurement patterns. 
For TQC the relation between the measurement and the desired gate is already transparent. 

Our first way of relating TQC and 1WQC was proposed by Aliferis and Leung JH]- The 
basic idea is to identify suitable pairs of consecutive 1-qubit measurements in 1WQC with a 
Bell measurement of a teleportation. This approach is strongly suggested by patterns such 
as the one in figure El With reference to the basic teleportation scheme in figure ^ we note 
the following. CX transforms Bell states into product states: 

cxiBij) = \(-iy)\i) 

(where the first ket is a |db) state according to the given sign). Hence the Bell measurement 
on 12 can be performed by the entangling operation CX followed by the 1-qubit measure- 
ments X\ and Z%. Similarly the |I?oo) state of 23 is CX |+) 2 |0) 3 which we can alternatively 
write as CX |+) 3 |0) 2 (where subscripts denote the qubit number and we adopt the nota- 
tion for the asymmetrical CX operation that the first qubit listed is the control qubit). 
Thus teleportation fully decomposes into entangling operations and 1-qubit measurements 
as shown in figure 1121 So, in view of figure El teleportation can be interpreted as: start 
with \ip)i |0) 2 |+)3, entangle suitably with CX's, then do X\ and Z2 measurements, which is 
structurally just like the 1WQC paradigm. But we have a slight mismatch in choice of primi- 
tives: TQC uses CX and |-Bqo) states whereas 1WQC is based on CZ and |_ff) = CZ |+) |+) 
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Figure 12: Teleportation as entangling operations and 1-qubit measurements. The diagram is read 
from left to right. The vertical lines denote CX operations with control and target marked as • and 
X respectively. Note that the two CX operations, with the shown opposite orientations, commute. 
Alternatively we could have put H on line 2 and the first CX the other way up, but then the two 
CX's would not commute. 



states, so the correspondence involves a sprinkling of Hadamard operations to interconvert 
these ingredients. 

Extending this idea, we find that other pairs of of measurements in 1WQC (such as 
X\M{d)2 and M\{9)X2) can be interpreted as rotated Bell measurements, but only special 
pairs of such consecutive 1-qubit measurements can be fused together to form Bell mea- 
surements. We refer to |15j for further details that we will not need here. Although we 
are able to reconstruct rotated Bell measurements for the 1WQC implementation of a full 
universal set of gates, this interpretation of 1WQC has the drawback that single 1-qubit 
measurements individually cannot be interpreted in terms of TQC. 

Our second relationship between TQC and 1WQC, proposed by Childs et al.^H] and 
Jorrand et al.|17j. is a further development of the ideas in figure IT2l As noted in the caption, 
the two CX operations commute. Also the X\ measurement commutes with all subsequent 
operations on 23. Thus we can change the order of actions to: (entangle 12, measure X\), 
then (entangle 23, measure Z2), obtaining a sequence of two operations of the same form 
viz. (entangle 12, measure 1) to obtain a state at 2. As noted in figure [71 it is exactly this 
kind of operation that drives 1WQC too, so we can regard it as a common fundamental 
primitive underlying both models (and sometimes called "one-bit teleportation"). 

A slight awkwardness in figure^] is the lack of uniformity of actions: the CX's act in 
opposite orientations (necessary for commutativity) and we have a single H gate as well. 
But this can be easily remedied: instead of the usual Bell state based teleportation scheme 
we consider teleportation with maximally entangled state \H) = CZ\+) |+) at 23 and its 
associated Bell measurement on 12 given by the basis 

{I ® I \H) ,X ® I \H) , Z ® I \H) ,XZ ® I \H)}. (3) 

Then note that CZ maps this basis to {|+) |+) , |+) |— ) , |— ) |— ) , |— ) |+)} so the Bell mea- 
surement is equivalent to applying CZ and then measuring X\ and Xi- Also unlike CX, 
CZ is symmetrical so the picture as in figure ED for this teleportation process is fully uni- 
form. The two "one-bit teleportations" are now identical, in fact corresponding exactly to 
the process in figure El implemented sequentially. 
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Figure 13: Valence bond solid or matrix product state for the cluster state. Each "valence bond" 
line denotes the maximally entangled state \H) of two qubits. Each site (circled) has 4 or 3 or 
2 qubits. At each site we consider the two dimensional subspace spanned by the two kets of "all 
zeroes" and "all ones" . 



7.1 Matrix product state relationship of TQC and 1WQC 

Our third and most remarkable connection between TQC and 1WQC, proposed by Ver- 
straete and CiracfPB]. is based on the formalism of so-called valence bond solids or matrix 
product states. In this correspondence each single 1-qubit measurement of 1WQC will be 
interpreted in terms of a full single teleportation. 

Consider a 2-dimensional grid of states \H) = CZ\+)\+) as shown in figure El Let 
| grid) denote the total state of all the qubits. At each site in the figure consider the 
two dimensional subspace spanned by {|00 . . . 0) , |11 . . . 1)} and the associated projector, 
renaming these two basis states as |6) and l): 

n= |6) (oo...o| + |i)(ii...i|. 

Applying II to (grid) we obtain a state with a single qubit at each site (and subnormalised 
because of the projection). 

Lemma 3: The multi-qubit state LT |grid) (after normalisation) is precisely the 1WQC 
cluster state. 

Proof: We first note the fact of figure ITU Now consider a 1-dimensional grid of \H) states 
as in figure IT51 Apply II at each node. Each bond is already CZ\+) |+) by definition. By 
the fact in figure El the projections at sites 23 and 45 simply serve to apply CZ between 
qubits 2 and 5. Hence the whole projected state is just CZ applied to all connecting pairs 
in |+) |+) ... |+) i.e. the 1-dimensional cluster state. This argument easily generalises to 
the 2-dimensional geometry of figure IT31 D 

Next consider using \H) states for TQC via application of rotated versions of the as- 
sociated basic Bell measurement eq. ©• For clarity of the essential idea, consider the 
1-dimensional case of figure IT5l Let us calculate the rotated Bell basis corresponding to the 
1-qubit gate 

1 / 1 e~ i6 \ 
^ = VS ( 1 ) ■ 
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Figure 14: Dots denote qubits and the connecting line denotes the maximally entangled state \H). 
If we project each circled site to span{|00) , |11)} then the resulting state (after normalisation) is 
GZ\ih)\ih). 



• • • • • » 

2 3 4 5 6 7 



Figure 15: One dimensional valence bond solid. 



The basis is given by 

\a) = {W(-6)^o a )®I\H) 

where a = 0, 1, 2, 3 and do = I, o\ = X, 02 = Z , = XZ . A direct calculation gives the 
first two states as 

|a = 0) = ^(|00)+e ie |ll» 
|o = l> = ^(|00>-e«|ll» 

and the remaining two span the orthogonal complement of span{|00) , |11)} at the site. 
Thus remarkably these first two Bell states lie within the 1-qubit subspace determined by 
the Il-projection. Furthermore this part of the Bell measurement corresponds precisely to 
the measurement basis M(9) = {|o) ±e l9 |l)} on the projected site i.e. on the cluster state. 
(Note that for other more general 1-qubit gates U, the corresponding rotated Bell basis 
states do not generally lie in a simple way relative to the Il-projected subspaces.) 

Stated otherwise, M(9) measurements on the cluster state (which is the basic ingredient 
of 1WQC, c.f. figure [7J) can be thought of as teleportations in the TQC formalism, where 
the teleportations always produce one of the a = 0, 1 outcomes (and not a = 2, 3) i.e. 
teleportations that have been "cut down" by the II projection. In a similar way all other 
1WQC ingredients (viz. Z measurements and the CZ measurement pattern) can be seen 
as descendants under the IT projection of teleportations on a valence bond grid state (grid). 
We omit further details which may be found in |18j . 
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8 Measurement based models and computational complexity 



The gate array model of quantum computation provides a transparent formalism for the 
theoretical study of quantum computation and its computational complexity features com- 
pared to classical computation. So why should we bother with further exotic models such 
as the measurement based models? Indeed our measurement based models are readily seen 
to be polynomial time equivalent to the gate array model i.e. each model can simulate the 
other with only a polynomial (i.e. modest) overhead of resources (number of qubits and 
computational steps). To see this first recall that the standard gate array model (allowing 
measurements only at the end and only in the Z basis) can be easily generalised to allow 
measurements along the way with subsequent choices of further gates and measurements 
being allowed to depend on earlier measurement outcomes. Indeed consider a measurement 
in a basis {U |0) , U |1)} on a qubit B applied during the course of a gate array process. To 
regain the standard gate array paradigm, for each such measurement we adjoin an extra 
ancillary qubit A, initially in state |0) and replace the measurement by the following: apply 
W to B and the apply CX to qubits BA. This simulates a coherent representation of the 
measurement in which qubit A plays the role of a pointer system. Subsequent gates that 
depend on the measurement outcome are replaced by a corresponding controlled operation, 
controlled by the state of A (written in the Z basis). In this way we purge all intermedi- 
ate measurements from the body of the array and a measurement of each ancilla in the Z 
basis at the end results in a standard gate array process which is equivalent to the given 
non-standard one. 

Using the above technique any 1WQC process is easily converted into an equivalent 
(standard) gate array process. We first build the required cluster state using an array 
of CZ gates acting on |+) = H\0) states and introduce an ancilla A for each 1-qubit 
1WQC measurement. For each M(9) measurement we introduce an extra gate W\6) which 
transforms the M(9) basis to the standard basis. The overhead in number of qubits and 
gates in this simulation is at most linear. 

Conversely given any gate array (based say on one of our previously considered universal 
sets of gates) we have seen how it can be translated into a measurement pattern on a suitably 
large cluster state. If K is the largest size of the 1WQC measurement pattern for any gate 
in our universal set then the number of qubits and computational steps increases by at most 
a factor of K i.e. the resource overhead is again linear. 

Polynomial time equivalence of computational models is important in computational 
complexity theory because such models have the same class of polynomial time computa- 
tions. But polynomial time equivalence does not preserve more subtle structural features 
of computations, such as parallelisability. Indeed already in the context of classical com- 
putation it is well known that the (one tape) turing machine model is polynomial time 
equivalent to the (classical) gate array model yet the turing machine model does not even 
have a natural notion of parallalisability at all, whereas the gate array model does! (i.e. 
doing gates simultaneously in parallel). 

In contrast to the quantum gate array model, the formalism of measurement based 
models offers new perspectives for parallelisability issues. We have already noted the fun- 
damental feature that measurements on different subsystems of an entangled state always 
commute so long as the choice of measurement is not adaptive i.e. not dependent on the 
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outcome of another measurement. We have seen examples of processes which are inherently 
sequential for gate arrays (e.g. sequences of Clifford gates) that become parallelisable in 
the measurement based models. 

The measurement based models have a further novel feature: they provide a natural 
formalism for separating a quantum algorithm into "classical parts and quantum parts". 
In contrast, in the gate array model every computational step is viewed as being quantum. 
The notion of classical-quantum separation becomes more compelling when we consider 
say, Shor's algorithm in its full totality, including the significant amount of non-trivial 
classical post-processing of measurement results needed to reach the final answer. It seems 
inappropriate to view this post-processing as a quantum process (albeit one that maintains 
the computational basis)! 

In measurement based computation the quantum parts of the algorithm are the quan- 
tum measurements done in parallel layers and the interspersed classical parts correspond 
to the adaptive choices of measurement bases, determined by classical computations on the 
previous layers' measurement outcomes. We may generalise this formalism in the following 
way: we allow (adaptively chosen) unitary gates as well as measurements within the quan- 
tum parts. We allow quantum layers to have only depth 1 (so a depth K quantum process 
is regarded as K layers with no interspersed classical computations) whereas classical layers 
can have any depth i.e. we are less concerned about controlling their structure. 

In this formalism any quantum computation is viewed as a sequence of classical and 
quantum layers. The total quantum state is passed from one quantum layer to the next and 
the quantum actions carried out in the next layer are determined by classical computations 
on measurement outcomes from previous layers. 

Any polynomial time quantum computation (say in the gate array model) can clearly 
be implemented with a polynomial number of quantum layers (and no interspersed classical 
layers) but the above formalism suggests a novel structural conjecture: 
Conjecture: Any polynomial time quantum algorithm can be implemented with only 
O(logn) quantum layers interspersed with polynomial time classical computations. □ 

This conjecture, asserting an exponential reduction in the essential "quantum content" 
of any quantum algorithm, has no analogue in classical complexity theory (where there 
is no notion of classical-quantum separation). Intuitively we are conjecturing that poly- 
nomial time classical computation needs relatively little "quantum assistance" to achieve 
the full power of polynomial time quantum computation. Although the conjecture remains 
unproven in general, we note that Cleve and Watrous JH] have shown that it holds true for 
Shor's algorithm. 
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